Anomalous charge transport in triplet superconductor junctions 
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Charge transport properties of a diffusive normal metal /triplet superconductor (DN/TS) junction 
are studied based on the Keldysh-Nambu quasiclassical Green's function formalism. Contrary to the 
unconventional singlet superconductor junction case, the mid gap Andreev resonant state (MARS) 
at the interface of the TS is shown to enhance the proximity effect in the DN. The total resistance 
of the DN/TS junction is drastically reduced and is completely independent of the resistance of 
the DN in the extreme case. Such anomalous transport accompanies a giant zero-bias peak in the 
conductance spectra and a zero-energy peak of the local density of states in the DN region. These 
striking features manifest the presence of novel proximity effect peculiar to triplet superconductor 
junctions. 



Physics of superconducting junctions has been one of 
the exciting fields of solid state physics in this decade . 
In diffusive normal metal / conventional singlet s-wave 
superconductor (DN/CSS) junctions, it is known that 
the phase coherence between an incoming electron and 
an Andreev reflected hole plays an essential role in caus- 
ing the proximity effect in the DN[1|1]. It is also known 
that the total resistance R of DN/CSS junctions does not 
follow the simple Ohm's rule, that is, R = Ro -\- Rr^^o 
where Rjj is the resistance in the DN and Rrd=o is the 
resistance of the DN/CSS interface. The reduction of R 
becomes prominent for low transparent junctions, where 
R < Rd -\- Rro=o is satisfied 2]. The lower limit of R 
is given by Rd + -Ro/2 using the Sharvin resistance i?o, 
where Rr^=o = Rpf^ is satisfied only for perfect trans- 
parent junctions Previous investigations of the prox- 
imity effect, however, are limited to DN/CSS junctions. 
Stimulated by the successive discovery of unconventional 
superconductors, we are tempted to expect novel proxim- 
ity effect in junctions with unconventional pairing, e.g., p- 
wave, and rf-wave, where pair potentials have sign change 
on the Fermi surface. 

In unconventional superconductor junctions, reflecting 
the internal phase of the pair potential, charge transport 
becomes essentially phase sensitive. The most dramatic 
effect is the appearance of zero bias conductance peak 
(ZBCP) 0, 01 in tunneling spectroscopy due to the for- 
mation of the mid gap Andreev resonant state (MARS) 
. The origin of the MARS is due to the anomalous in- 
terference effect of quasiparticles at the interface, where 
injected and reflected quasiparticles feel different sign of 
the pair potentials |^ . It is an interesting issue to clarify 
the role of the MARS on the transport properties of su- 
perconducting junctions. Recently, we have developed a 
theory of proximity effect, which is available for a diffu- 
sive normal metal / unconventional singlet superconduc- 
tor (DN/USS) junction, where USS indicates anisotropic 
singlet pairing like d-wave 0,13 ■ Unfortunately, however, 
it is revealed that the proximity effect and the MARS 



compete with each other in DN/USS junctions. Although 
the interface resistance Rr^=o is reduced by the MARS 
irrespective of the magnitude of the transparency at the 
interface, the resulting R is always larger than Ro/2+Rd. 
This is because the angular average of many channels at 
the DN/USS interface destruct the phase coherence of 
the MARS and the proximity effect (see Fig. 1). This 
destructive angular average is due to the sign change of 
the pair potentials felt by quasiparticles with injection 
angle (j) and those with —(p, where the angle (j) is mea- 
sured from the direction normal to the junction interface 
(see Fig. 1). However, in diffusive normal metal / triplet 
superconductor (DN/TS) junctions, we can escape from 
the above destructive average (see Fig. 1). We can ex- 
pect enhanced proximity effect by the MARS. In order to 
study this significantly novel charge transport, we must 
construct a novel theor y fo r DN/TS junctions beyond 
preexisting ones 0, H, 0,10 ■ This is in fact very timely 
since triplet superconductors have been discovered suc- 
cessively very recently 

In the present paper, we derive a conductance formula 
for DN/TS junctions based on the Keldysh-Nambu (KN) 
Green's function formalism 0I3- The total zero voltage 
resistance R in the DN/TS junctions is significantly re- 
duced by the enhanced proximity effect in the presence 
of the MARS. At the same time, local density of states 
(LDOS) in the DN region has zero energy peak (ZEP) 
due to the penetration of the MARS into the DN re- 
gion from the triplet superconductor (TS) side of the 
DN/TS interface. It is remarkable that when R^ is suffi- 
ciently larger than the Sharvin resistance Rq, R is given 
by i? = Rq/C-, which can become much smaller than 
the preexisting lower limit value of R, i.e., Ro/2 + Rjj. 
In the above, C_ is a constant completely independent 
of both Rd and Rb, where Rb denotes the interface re- 
sistance in the normal state. When all quasiparticles in- 
jected at the interface feel the MARS, R is reduced to 
be R = i?o/2 irrespective of the magnitude of Rd and 
Rb. The line shape of the bias voltage V dependent con- 
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ductance a{eV) has a giant ZBCP. These novel features 
have never been expected either in DN/CSS or DN/USS 
junctions. 

We consider a DN/TS junction with TS terminal 
and normal reservoir (N) connected by a quasi-one- 
dimensional diffusive conductor (DN) having a resistance 
Rd- The flat interface between the DN and the TS has 
a resistance Rb while the DN/N interface has zero re- 
sistance. The positions of the DN/N interface and the 
DN/TS interface are denoted as a; = —L and a; = 0, re- 
spectively as shown in Ref . Q . We restrict our attention 
to triplet superconductors with Sz — that preserves 
time reversal symmetry. Sz denotes the z component 
of the total spin of a Cooper pair. It is by no means 
easy to formulate a charge transport of DN/TS junc- 
tions since the quasiparticle Green's function has no an- 
gular dependence by the impurity scatteringin the DN. 
However, as shown in our previous paper 0, 1^, if we 
concentrate on the matrix currents ^3 '^i^ the TS 
to or from the DN, we can make a boundary condition 
of the KN Green's function. We assume that the con- 
striction area between the DN and the TS is subdivided 
into an anisotropic zone in the DN, two ballistic zones in 
the DN and the TS, and a scattering zone, where both 
ballistic and diffusive regimes can be covered 0. The 
sizes of the ballistic zone in the DN and the scattering 
zone in the current flow direction are much shorter than 
the coherence length 0; S 0- The scattering zone is 
modeled by an insulating delta function barrier with the 
transparency T{(j)) = 4cos^ 0/(4cos^ </> + Z'^), where Z 
is a dimensionless constant and (j) is measured from the 
interface normal to the junction 8]. The boundary con- 
dition for the KN Green's function in the DN [Gn{x)] at 
the DN/TS interface is given by. 
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{I): (1) 



using matrix current / |3, la| • Average over the angle of 
injected particles at the interface is defined by 

(/(^)) = / dcj) COS (j) I {(f) / d<pT{(j)) COS (j) (2) 

' J-tt/2 

with I{4>) = /. The resistance of the interface Rb is given 
by Rb =< Ro >■ As shown in the eq. (2) in our previous 
paper 0, matrix current / is a function of T{<p), Gi — 
Gn{x = 0_), G2+, and G2-, where G2+ {G2-) denotes 
the outgoing (incoming) Green's function in the TS. The 
Green's functions are fixed in the "TS" terminal and in 
the "N" terminal, and the voltage V is applied to the 
"N" terminal located a.t x = —L. Gn{x) is determined 
from the Usadel equation with eq. (1). If we denote 
the retarded part of GNix) and G2± as i?Ar(a:) and i?2± 
0, the following equations are satisfied, Rj^{—L) = f^, 
R2± = {f±fy+g±Tz) with /± = A±(</))/JA|(^)-e2 



and g± = e/y — A j_ (</)), using the Pauli matrices, e 
denotes the energy of the quasiparticles measured from 
the Fermi energy. A+(0) [A_ (</>)] is the pair potential 
felt by the outgoing (incoming) quasiparticles (see Fig. 
1). After some algebra, we can show that Rn{x) is given 
by sin9{x) ■ + cos9{x) ■ f^. The spatial dependence of 
d{x) in the DN is determined by the following equation 



D-^9{x) + 2iesm[e{x)] = 0, 



(3) 



with diffusion constant D in the DN. Taking the retarded 
part of Eq. we obtain 



L dB{x) 



Rd dx 



Rb ' 



(4) 



Br 



(riCos6>Q - Fa sm0o)T{(f>) 
(2 - r(0))F3 + T((^)[cos0or2 + smOaTi 



with 00 = 9{x = 0_), F2 = g++g-, T3 = l + f+f^+g+g_ 
and Fi = i{f+g- —g+f~)- Here, we focus on e = 0, where 
the left hand side of Eq. |0J is reduced to be Oq/Rd- 
We define F±((/)) as F±((/)) = lim,^o/± = sign(A±(0)) 
with sign(A+((/))) = 1(-1) for A+(0) > 0(< 0). In the 
following, we define the terminology as follows : when 
A A -{(f)) > is satisfied, quasiparticles are in the 
conventional channels (CC), while when A+(0)A_(0) < 
is satisfied, quasiparticles are in the unconventional 
channels (UC). In UC, quasiparticles feel the MARS 
while in CC quasiparticles do not feel the MARS. Follow- 
ing the above definition, F^{±<j)) = F_(±0) is satisfied 
for CC, while for UC, F+(±0) = -F_(±(/)) is satisfied. 
From cq. (4), we can show that is zero for CC and 
Bfj = iF+{(t)) for UC, respectively. Then, 6q becomes 



9q = iRoG-fRo, C_ 



F+{(j)) cos (j)d(j), (5) 



UC 



where jjj^ means the (f) integral only from the UC within 
— 7r/2 < (j> < 7r/2. A remarkable feature is that Oq be- 
comes a purely imaginary number as shown below. From 
eq. (3), 9{x) at e = becomes 6{x) = {x + L)6o/L. 
Since the LDOS of the quasiparticles in the DN region 
renormalized by its value in normal state is given by 
p(e) = Real[cos0(a:;)], p(0) is always larger than unity 
except for x = —L. This means that p(0) is enhanced 
due to the enhanced proximity effect by the MARS and 
has a ZEP as shown later. In the preexisting theories of 
DN/CSS and DN/USS junctions, 6*0 at e = is always a 
real number and p(x) never exceeds unity. 

In order to understand the essential difference between 
the DN/TS junctions and the DN/USS junctions intu- 
itively, we consider four simplified cases, (a)the DN/CSS 
junctions with the CC, (b)the DN/CSS junctions with 
the UC, (c)the DN/TS junctions with the CC, and (d)the 
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DN/TS junctions with the UC for all </>. This situation is 
actually realized by choosing d-wave pair potential with 
A±(0) = Ao cos[2((^ =F a)] (Figs. 1(a) and 1(b)) and 
p-wave pair potential with A±((f>) — ztAp cos(0 =F a) 
(Figs. 1(c) and 1(d)) as a prototype of USS and TS, 
respectively. The relations F_|_(±(/)) = F_(=F0) and 
F+{±(t>) = -F_{t<P), are satisfied for the DN/USS and 
the DN/TS junctions, respectively both for the UC and 
for the CC. For the DN/USS junctions with the CC, since 
F±{4') = F±{—4>) is satisfied, the contribution to is 
not cancelled by angular averaging [Fig. 1(a)]. For the 
DN/USS junctions with the UC and the DN/TS junc- 
tions with the CC, since F±{(f) = —F±{—(j)) is satisfied, 
the contribution to Oq is cancelled by the angular average 
[Figs. 1(b) and 1(c)]. However, for the DN/TS junctions 
with UC, since F±{(j)) = F±{—ip) is satisfied, the contri- 
bution to 6*0 is non-zero [Fig. 1(d)]. It is remarkable that 
for the DN/TS junctions the CC do not contribute to the 
proximity effect while UC do. 
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FIG. 1: Trajectories in the scattering process for outgoing 
(incoming) quasiparticles for the DN/USS and the DN/TS 
junctions and the corresponding A±{(j)) are schematically il- 
lustrated. As a prototype, we choose A± = Aocos[2(^ ^ 
a)] for the DN/USS junctions [(a) and (b)] and A± = 
±Ao cos[(0=pq)] for the DN/TS junctions [(c) and (d)], respec- 
tively, where Ao is the maximum value of the pair potential. 
The measure of the proximity effect is determined by the 
integration for all injection angles, i.e., the angular average 
over the hatched area. 

The total zero voltage resistance of the DN/TS junc- 
tion can be given as 



R 



Rd 

17 



dx 



Rt 



cos^ 



{x) {ho) ■ 



(6) 



where ho is obtained from the Keldysh part of / and is 
a complex function of g±, f±, 6o, and T{(j)). At e = 0, 
Ibo becomes 1 -I- exp(2 | 6*0 |) for UC, while for CC, Ibo 



becomes ^— 

by 



[2-T(0)]3 '"'^^^^o- After simple algebra, R is given 



R — Ro{ 



tanh6'o'i 

C- 



[l + exp(2|0oz|)]C++coshX,D 



} 



with = - RdC^/Ro, C+ = /, 



(7) 



UC 



cos 



and 

D — Jij^ [2-T(^)\i '^^^ where J^^ means the (j) integral 
only from the CC. The resulting R is given by Rq/C^ for 
sufficiently large Rd/Ro and is independent both of Rd 
and Rb except for the very special case with C_ = 0. The 
magnitude of R can become much smaller than the preex- 
isting lower limit value of i?, i.e., i?o/2 + Rd- This giant 
reduction of R is due to the enhanced proximity effect by 
the MARS. When all quasiparticles feel the MARS inde- 
pendent of (j), C_ is given by 2 and R becomes R = Rq/2 
independent of Rd and Rb. This interesting situation 
is actually realized for a Pa;-wave case, where pair po- 
tentials are given by A±((/)) — ±Aocos(/) [see case Fig. 
1(d)]. Above enhanced proximity effect by the MARS 
is a completely novel proximity effect, which has never 
been expected in any preexisting theories 0, S S IS 0] ■ 

In order to understand this novel proximity effect in 
detail, we focus on the LDOS p(e) in the DN region nor- 
malized by its value in the normal state. We choose Px- 
wave pairing as a prototype of the TS. As a reference, 
we compare the results with those for the DN/USS junc- 
tion with c?2;y-wave pair potential, where lS.±{<j)) is given 
by A±(^) = ±Aosin(20). Although p(e) at the TS side 
of the DN/TS interface and that at the USS side of the 
DN/USS interface both have ZEP by the formation of 
the MARS, p(e) in the DN has a drastic difference be- 
tween the two cases. For the p^-wave case, the LDOS 
in the DN region has a ZEP. The height of the ZEP 
p(0) is given by cosh[2i?£)(a; -|- L)/{LRq)\ and the order 
of its width is Exh- On the other hand, for the dxy- 
wave case, p(e) is always unity independent of Rd due 
to the absence of the proximity effect. Contrary to the 
DN/TS junction, the MARS formed at the USS side of 
the DN/USS interface can not penetrate into the DN. As 
seen from p(e) in DN/TS junctions, the significant reduc- 
tion of R originates from the penetration of the MARS 
into the DN. Although we have actually shown the ex- 
istence of ZEP of p(e) in the DN for px-wa.ve case as a 
prototype, ZEP is universally expected for DN/TS junc- 
tions with the MARS at the interface independent of the 
detailed shape of the pair potentials of the TS. On the 
other hand, the LDOS in the DN region of DN/USS junc- 
tions do not have ZEP. Using this clear and qualitative 
difference of LDOS in the DN region between the DN/TS 
and DN/USS junctions, we can identify the triplet sym- 
metry of the pair potentials. 

Finally, we focus on the line shape of the tunneling 
conductance of DN/TS junctions for non zero voltage, 
which is defined by a{eV) = RB/R{eV). We choose 
Paj-wave pairing as a prototype. As a reference, we com- 
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FIG. 2: p(e) is plotted for the DN/TS junction with p^-wave 
superconductor (left panel) and the DN/USS junction with 
dxy-wave superconductor (right panel) for Z = 1.5, Rd/Rb ~ 
0.5, and Eth = 0.02Ao with various x in the DN. a; x = 0_, 
b: x = — 1//4 and c: x = —L/2. 

pare the results with those for the DN/USS junction with 
dajy-wave pair potential. Although both p^^-wave and 
da;j,-wave cases have similar line shapes of the voltage- 
dependent conductance with the ZBCP as a function of 
eV for ballistic junctions ^^li Rd = case, we can 
classify these two for Rd / as shown in Fig. |3| We 
choose the Thouless energy Etii as Eth — 0.02Ao. For 
the pj;-wave case, since all injected quasiparticles feel the 
MARS, cr(0) = 2Rb/Ro is satisfied for any Rd- For 
Rd 7^ 0, a{eV) can be expressed by the summation of the 
broad ZBCP and the narrow one, where (Tn = Rq/Rb 
denotes the angular averaged transparency of the junc- 
tion. The width of the former one is proportional to 
cttvAq and the latter one is the Thouless energy. How- 
ever, with the increase of R^i/Rb, (j{eV) for | eV |> Exh 
is suppressed, and the ratio of ct(0) to its background 
value is largely enhanced. We can call this largely en- 
hanced a{Q) as giant ZBCP (curve & or c in the left 
panel of Fig. 3). By contrast, for dxy-wave case, since 
cr(0) ~ 2Rb / {Ro + IRd) is satisfied, (t(0) is reduced with 
the increase of Rd- The width of ZBCP is proportional 
to AoCTAr and is not changed with the increase oi Rd / Rb 
due to the absence of the proximity effect. 
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FIG. 3: a(eV) is plotted as a function of eV for the 
DN/TS junctions with pa;-wave superconductor (left panel) 
and the DN/USS junctions with dxy-wsye superconductor 
(right panel) for Z = 1.5 and Eth = 0.02Ao. a: Rd/Rb = 0, 
b: Rd/Rb = 0.5 and c: Rd/Rb = 1. 



In conclusion, we have presented a theory of charge 
transport in the DN/TS junctions. The total resistance 
R can become much smaller than the preexisting lower 
limit value of R for DN/CSS and DN/USS junctions, 
i.e., Ro/2 + Rd- As the extreme case where all injected 
electrons feel the MARS, R is reduced to be Ro/2. The 
significant reduction of R is due to the enhanced prox- 
imity effect by the MARS and the resulting LDOS in 
the DN region has a ZEP. At the same time, we can 
expect giant ZBCP. The above enhanced proximity ef- 
fect is a novel proximity effect, which has never been 
expected in DN/CSS and DN/USS junctions. We have 
shown in the present paper that these effects are actually 
realizable for DN/TS junctions with p-wave pair poten- 
tial. We believe that these features are easily verifiable 
in experiments since a mesoscopic interference effect due 
to the proximity effect has recently been observed in high 
Tc cuprate junctions Similar experiments are tech- 
nologically possible for junctions composed of Sr2Ru04, 
where triplet pairing is believed to be realized jllj] . 

The authors appreciate useful and fruitful discussions 
with Y. V. Nazarov, A. A. Golubov J. Inoue, Y. Asano 
and K. Kuroki. 
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